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We use non-perturbative linked-cluster expansions to determine the ground-state energy per site 
of the spin-one Heisenberg model on the kagome lattice. To this end, a parameter is introduced 
allowing to interpolate between a fully trimerized state and the isotropic model. The gronnd-state 
energy per site of the full graph decomposition up to graphs of six triangles (18 spins) displays a 
complex behaviour as a function of this parameter close to the isotropic model which we attribute 
to divergencies of partial series in the graph expansion of quasi-ld unfrustrated chain graphs. More 
concretely, these divergencies can be traced back to a quantum critical point of the one-dimensional 
unfrustrated chain of coupled triangles. Interestingly, the reorganization of the non-perturbative 
linked-cluster expansion in terms of clusters with enhanced symmetry yields a ground-state energy 
per site of the isotropic two-dimensional model being in quantitative agreement with other numerical 
approaches in favor of a spontaneous trimerization of the system. Our findings are of general 
importance for any non-perturbative linked-cluster expansion on geometrically frustrated systems. 
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I. INTRODUCTION 

The search for exotic phases of correlated quantum 
matter is one of the most attractive issues in physics, 
since one expects fascinating types of collective quantum 
behaviour hosting unconventional physical properties of 
potential interest for future technologies. To this end, the 
interplay of frustrated interactions and strong quantum 
fluctuations is believed to represent the key knobes one 
can tune to stabilize unconventional types of quantum 
order. 

In this context the physics of the quantum spin- 
1/2 Heisenberg model on the highly frustrated two- 
dimensional kagome lattice plays a paradigmatic role 
and a great amount of research has been performed over 
the last decades. The long history of that problem di¬ 
rectly reflects the enormous complexity of such highly 
frustrated q uant um magnets. Due to recent numerical 
investigation^U^I^ the nature of the ground state is most 
likely a quantum spin liquid exhibiting exotic topological 
order. Elementary excitations above the topologically- 
ordered ground state are therefore neither bosons are 
fermions, but anyons carrying fractional spin quantum 
numbers. 

The physics of the spin-one cousin of the Heisen¬ 
berg model on the kagome lattice has been focused on 
recentlj®2l_ Interestingly, several works point towards 
a spontanenous trimerization of the systeirPl^, i.e. the 
symmetry between up- and down-triangles is broken 
spontaneously and a singlet ground state with no long- 
range magnetic order is found. As a consequence of 
breaking this discrete symmetry, the spin-one Heisenberg 
model is expected to be gapped. This is indeed in agree¬ 
ment with recent density matrix renormalization groim 
(DMRG) calculations pointing towards a hnite spin gapP. 
Nevertheless, there are also competing numerical works 


which argue in fa vor of different ordered or disordered 
ground state^^*^, and it is therefore important to further 
clarify the physical properties of the spin-one Heisenberg 
model on the kagome lattice. On the experimental side, 
there are also several materials like m-MPYNN-BFii^ti^ 
or Ni 3 V 20 ^^ which seem to correspond microscopically 
to spin-one systems on the kagome lattice. 

If the spin-one Heisenberg model on the kagome lat¬ 
tice displays spontaneous trimerization, then it should 
be possible to deform the ground state of the Heisen¬ 
berg model adiabatically up to the fully trimerized limit 
of isolated triangles without encountering any quantum 
phase transition. In this work we check this scenario by 
calculating the ground-state energy per site in the ther¬ 
modynamic limit using non-perturbative linked-cluster 
expansions (NLCEs)p2®l, yye set up a full graph decom¬ 
position in terms of up-triangles of the kagome lattice and 
we perform exact diagonalization (ED) on these graphs. 
We hnd a surprisingly complex convergence behaviour of 
the ground-state energy per site. This is shown to orig¬ 
inate from a quantum critical point present in the one¬ 
dimensional subsystem of coupled triangles. The traces 
of this quantum critical point can be consistently avoided 
by reorganizing the NLCE of the two-dimensional prob¬ 
lem in terms of more symmetric clusters. This reorga¬ 
nized NLCE displays a strongly improved convergence 
and our results are fully consistent with a trimerized 
ground state of the spin-one Heisenberg model on the 
kagome lattice. On the technical side, the appearance of 
divergent partial series in non-perturbative linked-cluster 
expansions is expected to be a generic feature in geomet¬ 
rically frustrated systems. 

The outline of this paper is as follows. In Sect. [H] we 
introduce the trimerized version of the spin-one Heisen¬ 
berg model on the kagome lattice. Next, we give technical 
details on the non-perturbative linked-cluster expansions 
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FIG. 1. (Color online) a) A schematic illustration of the spin- 
one Heisenberg model on the kagome lattice. Spins are lo¬ 
cated on the vertices of the lattice highlighted as circles. The 
Heisenberg exchange on up- (down-) triangles is J (AJ) and 
illustrated as black (gray/cyan) lines. The black boxes mark 
the center of up-triangles, b) The effective triangular lattice 
of up-triangles is shown. Note that up- and down-pointing 
triangles (light/green and dark/red areas) in the effective lat¬ 
tice are topologically distinct, since they correspond to differ¬ 
ent plaquettes (hexagons and down-triangles) in the original 
kagome lattice. 

performed in this work in Sect. |III| and we present all 
physical results in Sect. |IV[ Finally, the most important 
outcomes of our work are briefly concluded in Sect. [V| 


II. MODEL 

Let us consider a trimerized version of the spin-one 
Heisenberg model on the kagome lattice as illustrated in 
Fig. [3 The microscopic Hamiltonian is given by 

n = j s,-Sj + Aj ’ (1) 

(hi)6V 

so that the first (second) sum runs over all nearest neigh¬ 
bor sites on up-triangles (down-triangles). The real pa¬ 
rameter A G [0,1] allows an interpolation between the 
limit of isolated up-triangles A = 0 and the isotropic 
Heisenberg model at A = 1. In the following we focus 
on an antiferromagnetic exchange constant J > 0 and 
we set J = 1. 

For A = 0, the fully trimerized system has a unique 
singlet ground state which is given as the product state 


of singlets on individual up-triangles. Indeed, in con¬ 
trast to its spin-1/2 cousin where an isolated antiferro¬ 
magnetic triangle is four-fold degenerate, the addition of 
three spin-ones allows to form a unique singlet which has 
the lowest energy. All excited states at A = 0 are gapped 
and can be understood via the remaining 26 states of an 
isolated triangle and their magnetic quantum numbers. 

If the isotropic model at A = 1 displays a spontaneous 
trimerization, then one expects that ground states at 
A = 0 and A = 1 are adiabatically connected and, conse¬ 
quently, no energy gap should close in the whole interval 
A G [0,1]. This scenario we explore in the following by 
applying NLCEs which we explain next. 


III. NLCE 

The essential idea of NLCEs is to exploit the linked- 
cluster theorem so that actual numerical calculations are 
done on finite linked clusters, but the final results are 
valid directly in the thermodynamic limit. Generically, 
NLCEs consists of three steps: i) choosing and generat¬ 
ing the families of clusters or topologically distinct graphs 
used in the linked-cluster expansion, ii) performing nu¬ 
merical calculations on graphs extracting the physical 
quantities of interest, and iii) determing the reduced con¬ 
tributions specific to each graph and embed these contri¬ 
butions into the infinite lattice. 

The details of the choice of graphs (i) will be given 
below. Concerning ii), we are using ED with the Lanc- 
zos algorithm to determine the ground-state energy Eg’' 
on each graph and for each value of A. Apart from the 
exponential increase of graphs with the number of sites 
iV, the memory needed for ED is the limiting factor of 
NLCEs. Practically, we consider graphs with at most 18 
spin-ones in our calculations. 

In the third step iii), the ground-state energy per site 
Co in the thermodynamic limit is expressed as 

eo = XI ^ (2) 

j 

where the sum runs over all linked graphs Qj. The integer 
number is the so-called embedding factor specifying 
the combinatorical number how often graph Qj can be 
embedded into the infinite lattice. The reduced contri¬ 
bution Cq specific to graph Qj results from subtracting 
all contributions of subgraphs in order to avoid double 
counting. One has 

Ared=A^- E • (3) 

In the following we detail the two different types of 
NLCEs we performed for the spin-one Heisenberg model 
on the kagome lattice. 
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FIG. 2. (Color online) Illustration of all graphs used in the 
full graph expansion up to Nn = 3 triangles. The left panel 
corresponds to the original kagome lattice while the right 
panel shows the same graphs in the effective lattice where up- 
triangles (shown in black) are replaced by supersites (black 
boxes). 


A. Full-graph decomposition in triangles 

We are aiming at describing the trimerized ground 
state for A S [0,1]. Is is therefore reasonable to set up a 
full-graph expansion not in terms of individual spin sites, 
but in terms of up-triangles (or equivalently of down- 
triangles) as effective supersites. The effective lattice of 
these supersites corresponds to a triangular lattice as il¬ 
lustrated in the lower panel of Fig. The up- and down- 
triangles of the effective triangles are topologically dis¬ 
tinct, since they correspond to hexagons and triangles in 
the original kagome lattice. Respecting this property, we 
performed a full graph expansion in the effective lattice 
for all graphs up to Wr = 6 supersites, i.e. 6 up-triangles 
in the kagome lattice consisting of 18 spins one. 

All graphs up to Ntr = 3 are illustrated in Fig. in 
the original kagome lattice (left panel) as well as in the 


effective triangular lattice (right panel). In total, there 
are 1, 1, 4, 9, 46, 228 topologically distinct graphs for Ntr 
from 1 to 6. Finally, we have determined the ground-state 
energy per site taking into account all graphs up 

to Ntr S 6}. 

B. Reorganized graph expansion 

Next we introduce a reorganized graph expansion 
which turns out to be well suited for the spin-one Kagome 
Heisenberg model as detailed below. The essential idea 
is to restrict the calculation to only a small subset of 
clusters which, however, are highly symmetric. 

In terms of the effective triangular lattice as shown 
in Fig. [^, the important building blocks of this lattice 
are the two inequivalent triangles as displayed in the up¬ 
per panel of Fig. We use these two objects as the 
central blocks for the reorganized graph expansion. In 
the simplest approximation we restrict the NLCE 
calculation to these two clusters plus the two-site chain 
graph which needs to be taken into account in order to 
avoid double counting. Note that one has to consider the 
single supersite graph which therefore corresponds for¬ 
mally to This scheme is then extended by adding 

only clusters made (symmetrically) of effective up- and 
down-triangles so that neighboring triangles always share 
edges and not corners as shown in Fig. In the follow¬ 
ing we numerate this reorganized graph expansion by the 
total number of triangles involved so that we determined 
with n G {0,1, 2,3,4}. Finally, one can further 
reorganize the graph expansion to even more symmetric 
clusters made of effective plaquettes where a plaquette 
consists of one effective down- and up-triangle. This ex¬ 
pansion is performed exactly when restricting to even n. 


C. Extrapolation 

Within the NLCE we get a finite series for the ground- 
state energy per site eo in the thermodynamic limit. For 
the limit Wr oo one expects to obtain the exact re¬ 
sult. In this context the eg^*'^ corresponds to the partial 
sum of all orders up to Wr- To increase the convergence 
of the series there is a powerful technique called “ series 
acceleration”. For our problem the Wynn algorithrrfiSllII 
is of particular interest. It takes the form 

= 0 (4) 

(5) 

-I- _ ----- (6) 

^n+1 

The parameter n runs from 1 to Wr in the initial step 
and is lowered by 1 in each Wynn step. In practize. 
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FIG. 3. (Color online) Illustration of all graphs used in the 
reorganized expansion except the one consisting just of a sin¬ 
gle supersite. The idea is to use the two inequivalent triangles 
of the effective triangular lattice as the building blocks of a 
NLCE calculation. 



FIG. 4. (Color online) Illustration of the one-dimensional 
spin-one triangle chain. Spins are located on the vertices of 
the lattice highlighted as black circles. The Heisenberg ex¬ 
change on up-triangles is J (black lines) while the unfrus¬ 
trated Heisenberg exchange is J between triangles forming 
dimers (shown as gray/cyan lines). In the limit J —> 0 one 
can derive an effective spin-one Heisenberg chain in second- 
order perturbation theory in J/J for the spin ones building 
the tips of the black triangles (dark/green circles). The effec¬ 
tive Heisenberg exchange Jefr is antiferromagnetic and illus¬ 
trated as dotted lines. 

one can also restrict to certain orders, e.g. one could only 
take into account even or one omits lowest-order con¬ 
tributions by setting For odd k the entries 

usually diverge, whereas they converge for even k. At the 
end the sequence e^_ 2 ; (for I = 0,1, 2,...) converges to 
eo- 

IV. RESULTS 

We start the discussion with a one-dimensional trian¬ 
gle chain which represents a one-dimensional subsystem 



FIG. 5. (Golor online) Spin gap A as a function of A for the 
one-dimensional triangle chain with periodic boundary con¬ 
ditions, J = 1, and number of triangles Ntv The dashed 
vertical line marks the location of the quantum critical point 
in the thermodynamic limit. Inset left: Minimal spin gap Ac 
as a function of l/Ntr- The dashed line represents a linear 
fit through the data points signaling a quantum phase tran¬ 
sition in the thermodynamic limit. Inset middle: Ratio Ac 
corresponding to the minimal value of the spin gap Ac as 
a function of l/Wr. The dashed line represents a linear fit 
through the data points signaling the location of the quan¬ 
tum critical point in the thermodynamic limit. Inset right: 
Second derivative df eo as a function of A. 

of the trimerized Kagome lattice. This model is essen¬ 
tially unfrustrated (in terms of the effective triangles) 
and, in addition, the full graph expansion becomes very 
simple due to the one-dimensional nature. This system 
is shown to possess a quantum critical point which has 
severe consequences for the full-graph expansion of the 
two-dimensional trimerized Kagome lattice presented af¬ 
terwards. Finally, we discuss the physical properties of 
the reorganized NLCE which displays a very good con¬ 
vergence behaviour. 

A. One-dimensional triangle chain 

We consider the one-dimensional spin-one triangle 
chain illustrated in Fig. The Hamiltonian is given by 

H = J ^ S,-S.-bJ ^ s,-S, , (7) 

{i,j)^A {i,j) Gdimer 

so that the first sum runs over all nearest neighbor sites 
on triangles while the second sum contains all nearest 
neighbor sites between adjacent triangles forming spin- 
one dimers. In the limit J = AJ —>■ 0, one is close to the 
trimerized limit as for the two-dimensional kagome coun¬ 
terpart, while in the limit J —>■ 0 with J finite the system 
realizes singlets on dimers between triangles while there 
are an extensive number of isolated spin ones building the 
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FIG. 6. (Color online) Ground-state energy per site of 

the one-dimensional triangle chain as a function of A calcu¬ 
lated by NLCE using a full graph expansion for different Nti 
and setting J — 1. Upper inset: Reduced contributions 
as a function of A for different A^tr- Lower inset: Reduced 
contributions | as a function of l/A'tr for fixed values of 

A in a double logarithmic plot. Dashed black line represents 
an algebraic fit for the displayed data with A = 0.8. 



FIG. 7. (Color online) Comparison of the ground-state energy 
per site eo as a function of A for J = 1 calculated by NLCE 
using a full graph expansion for different Ntr (shown as sym¬ 
bols) and by iPEPS (shown as solid black line) taken from 
Ref. E The gray/cyan diamonds correspond to the Wynn 
extrapolation where all contributions from Wr = 2 to 
Ntr = 6 are taken into account. 


tips of the up-triangles. If one performs second-order de¬ 
generate perturbation theory in this extensive subspace 
of isolated spins, one obtains up to an irrelevant con¬ 
stant the Hamiltonian of an effective antiferromag¬ 
netic spin-one Heisenberg chain 

Heff = + Jeff ^ S., • S, , (8) 

(hi) 

where Jeff = This systems is known to real¬ 

ize a symmetry-protected top ologica lly-ordered ground 
state with a finite Haldane gapP^lMl One therefore ex¬ 
pects at least one phase transition when varying the ratio 

J/J = A. 

We have used ED on finite triangle chains with peri¬ 
odic boundary counditions consisting of up to fVtr = 6 
triangles to study the ground-state phase diagram. We 
calculated the ground-state energy per site cq and the 
spin gap A as the difference between the ground states 
in the = 1 and = 0 sector. The corresponding 
data are shown as a function of A in Fig. [^setting J = 1. 
We find clear evidence for a breakdown of the trimerized 
phase around Ac « 0.8: i) the spin-gap scales to zero for 
this ratio and ii) the second derivative develops a 
resonance in the same A-regime fully consistent with a 
divergence in the thermodynamic limit. Interestingly, we 
do not observe any further signatures for phase transi¬ 
tions in the range A > 1 which implies that the ground- 
state phase diagram of the triangle chain consists of two 
distinct quantum phases separated by a quantum phase 
transition at Ac < 1. 

Next we turn to our NLCE results for the triangle 
chain. The NLCE recipe for one-dimensional systems be¬ 


comes particularly simple. Graphs are linked chain seg¬ 
ments containing up-triangles, so for each there 
is just one graph, and one has for all embedding factors 
ytjj — Furthermore, the NLCE ground-state energy 
per site in the thermodynamic limit can be expressed 
analytically as 


since all other subcluster contributions cancel out ex¬ 
actly. The corresponding results for up to Ntr = 6 

are displayed in Fig. One observes a smooth behaviour 
and a very good convergence for increasing values of 
iVtr. This is even better seen in the reduced contribu¬ 
tions shown in the upper inset 

of Fig.|^ They get smaller with increasing and larger 
with increasing A. No signatures from the quantum crit¬ 
ical point around Ac can be detected in these quantities 
(the same is true for 9feg which is not shown). This 
is very likely due to the fact that quantum fluctuations 
from the trimerized phase are optimally taken into ac¬ 
count with the chosen clusters while corresponding fluc¬ 
tuations of the phase present for A > Ac are only captured 
deficiently. 

Interestingly, there is an alternative way to grasp the 
quantum critical behaviour in the presented NLCE data 
for the triangle chain. In gapped phases one expects an 
exponential convergence of the NLCE data with increas¬ 
ing TVtr due to the finite correlation length which should 
be the case for all A except Ac. In contrast, at Ac (or 
close to Ac for any finite Atr) there should be an alge¬ 
braic convergence due to the diverging correlation length 
when enlarging Atr (see also Ref. [1^ . In the lower inset 
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of Fig. [^the corresponding data of the reduced contri¬ 
butions IcQ^edl displayed as a function of 1/iVtr for 
several fixed values of A in a double logarithmic plot. 
Clearly, only the data set for A = 0.8 « Ac is well de¬ 
scribed by an algebraic convergence and there are visible 
deviations for values of A being smaller or larger than Ac. 

To summarize this part, we have found that the one¬ 
dimensional spin-one triangle chain possesses a quantum 
critical point around Ac ~ 0.8, i.e. the trimerized phase 
is not stable up to the isotropic value A = 1. As a con¬ 
sequence of this quantum critical behaviour, the reduced 
NLCE contributions of the ground-state energy per site 
only decay algebraically with cluster size close to Ac. 

B. Kagome lattice 

We start by discussing our NLCE results for the 
ground-state energy per site obtained from a full 
graph expansion in terms of triangles for the spin-one 
Heisenberg model on the kagome lattice. The corre¬ 
sponding data is shown together with infinite Projected 
Entangled Pair States (iPEPS) results taken from Ref. [8] 
in Fig. Compared to the case of the triangle chain 
discussed in the last section, the convergence behaviour 
is clearly not as good with increasing Ntr- Curves with 
even and odd A^tr belong to two different families of data 
sets, since all curves with even (odd) Atr are above (be¬ 
low) the iPEPS data. The NLCE curves converge very 
slowly for larger values of A and there are visible anoma¬ 
lies as a function of A already before the isotropic point 
A = 1. A certain improvement is achieved by applying 
the Wynn algorithm. In order to include the contribu¬ 
tion of the highest order Atr = 6, one has to exclude the 
lowest-order contribution as explained in Subsect. |III C| 
Therefore, the last step of the Wynn extrapolation yields 
which is also shown in Fig. H 

This poor convergence is a consequence of the quan¬ 
tum critical behaviour of the unfrustrated triangle chain 
deduced in the last subsection. Indeed, all the trian¬ 
gle chain clusters of the one-dimensional triangle chain 
are also present in the full graph expansion in terms of 
triangles on the kagome lattice. These clusters are es¬ 
sentially unfrustrated and we therefore find that reduced 
contributions are negative. Physically, this means that 
these clusters lower the energy per site of the trimerized 
ground state, while the contributions of other clusters 
typically increase the energy due to the geometric frus¬ 
tration. We stressed in Subsect. IIV Al that reduced con¬ 
tributions of chain graphs decay only algebraically 

with Atr close to Ac « 0.8. For the two-dimensional 
kagome lattice, this becomes problematic due to the fact 
that the embedding factor of these graphs grows expo¬ 
nentially with Atr- Consequently, the product 
which enters the NLCE expression for the ground-state 
energy per site Eq. diverges when enlarging the cluster 
size of the chain graphs (see also Fig. for A = 1). This 
divergent partial series of the unfrustrated chain graphs 



FIG. 8. (Color online) Total contribution ^n- 

frustrated chain graphs as a function of Atr for A = 1 where 
1 /^*' is the embedding factor and the reduced contribu¬ 

tion to the energy per site. Lower inset Embedding factor 
of unfrustrated chain graphs as a function of Atr for 
A = 1. 



FIG. 9. (Color online) Comparison of the ground-state energy 
per site eo for J = 1 as a function of A calculated by NLCE 
using reorganized graphs of different size (shown as symbols) 
and by iPEPS (shown as solid black line) taken from Ref. [H 
Diamonds correspond to two different Wynn extrapolations. 
The open black (filled gray/cyan) diamonds belong to 
(®o^even) where all (even order) contributions are taken into 
account. Inset. Zoom of the same data close to A = 1. 


has to be compensated by the overall contribution of the 
majority of remaining frustrated graphs to get a finite 
total result for the ground-state energy per site. There¬ 
fore, the poor convergence seen in Fig.j^for the full graph 
expansion in terms of triangles should not come as a sur¬ 
prise when truncating the NLCE series by taking a finite 
Ntr- 

One obvious idea is to reorganize the graph expansion 
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in such a way that there are no divergencies in partial 
series resulting from quantum critical behaviour of one¬ 
dimensional unfrustrated subsystems. And indeed, us¬ 
ing the two inequivalent effective triangles as frustrated 
and more symmetric building blocks, as introduced in 
Subsect. |IIIB[ the convergence of the obtained NLCE 
ground-state energy per site is dramatically improved as 
shown in Fig. Again, the results alternate for 

n being odd or even, but all curves smoothly approach 
the iPEPS result and absolute differences between NLCE 
and iPEPS are considerably smaller in the whole interval 


A S [0,1]. 

Further improvement is achieved by applying the 
Wynn algorithm. In order to include the contribution 
of the highest order Ntr = 4, we include the contribution 
of a single supersite as the lowest-order contribution (see 
also Subsect. |IIIC|. The last step of the Wynn extrap- 

For this extrapolation we find 


= ( 4 ) 


olation then yields Cq 
that the difference to the iPEPS data is similar com¬ 
pared to the one with the bare NLCE result at 

A = 1 (but with a different sign). Interestingly, if we per¬ 
form the Wynn algorithm restricting to even-order con¬ 
tributions, which corresponds to an effective plaquette 

( 2 ) 

expansion, the obtained extrapolant Cq shows almost 
perfect agreement with the existing numerical data even 
at A = I: The value of the extrapolated iPEPS ground- 
state energy per site is —1.4116(4) J from ReflHl the one 
calculated by DMRG from Ref. His —1.410(2) J while we 
find e'^Ln = -1.4114J. 


V. CONCLUSIONS 

In this work we have studied the trimerized spin-one 
Heisenberg model on the kagome lattice applying non- 
perturbative linked-cluster expansions. To this end a pa¬ 
rameter A is introduced which interpolates between the 
limit of isolated triangles and the isotropic Heisenberg 
model. We used a NLCE in terms of triangles for the 
ground-state energy per site which, by construction, be¬ 
comes perturbatively controlled for small A. 

This full graph expansion shows a non-monotonic 
convergence behaviour when approaching the isotropic 
Heisenberg model. Our findings can be traced back to 
quantum critical behaviour from unfrustrated quasi-Id 
chain graphs. For these graphs the embedding factor di¬ 


verges faster (exponentially) with increasing graph size 
than the suppression of reduced contributions (algebraic 
behaviour) leading to an overall divergence of this partial 
subclass of graphs. Such divergencies have to be compen¬ 
sated by the remaining majority of frustrated graphs in 
order to yield a physically meaningful result which is a 
subtle issue when truncating the sequence in a NLCE 
calculation. The appearance of divergencies in partial 
series of unfrustrated graphs is likely a generic feature 
of NLCE calculations on frustrated systems using a full 
graph expansion. 

In our opinion there are two possibilities to overcome 
this finding. First, one can reorganize the graph expanion 
so that it is better suited for the problem under study. 
More concretely, one must avoid divergent partial series 
of graphs. For the kagome spin-one Heisenberg model we 
have shown that a reorganization of graphs in terms of 
effective triangles (hexagons and triangles in the original 
kagome lattice) leads to an enormous improvement in 
the convergence giving quantitative agreement with other 
recent approaches. It would also be interesting in future 
to investigate other reorganized graph expansions for this 
problem. 

Second, one might wonder whether there exist an 
optimized full graph expansion which goes beyond the 
paradigm of using ED on graphs as recently deduced for 
NLCEs when targetting excited state^^. It was shown 
in Ref. [50] that graph-based continuous unitary transfor¬ 
mations (gCUTs) are capable to acount for the reduced 
symmetry of graphs which leads to artificial processes 
on graphs which are not present in the thermodynamic 
limit. It would be fascinating if a similar progress can be 
found for ground-state properties of frustrated systems. 

Finally, it would be interesting to determine the prop¬ 
erties of excited states in the spin-one kagome Heisenberg 
model which we leave for future studies. 
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